Characterizing multipartite symmetric Dicke states under the effects of noise 



On 
O 
O 
<N 

a 
a*. 



> 

m 
cn 

o 



S. Campbell, M. S. Tame, and M. Paternostro 

School of Mathematics and Physics, Queen's University, Belfast BT7 INN, United Kingdom 

(Dated: June 26, 2009) 

We study genuine multipartite entanglement (GME) in a system of n qubits prepared in symmetric Dicke 
states and subjected to the influences of noise. We provide general, setup-independent expressions for experi- 
mentally favorable tools such as fidelity- and collective spin-based entanglement witnesses, as well as entangled- 
class discriminators and multi-point correlation functions. Besides highlighting the effects of the environment 
on large qubit registers, we also discuss strategies for the robust detection of GME. Our work provides tech- 
niques and results for the experimental communities interested in investigating and characterizing multipartite 
entangled states by introducing realistic milestones for setup design and associated predictions. 
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The ability to classify entangled states Q, Q, H] and 
quantify their degree of correlation [4] is advancing to- 
gether with the capacity to experimentally produce interest- 
ing and useful forms of quantum correlated many-body sys- 
tems. The range of experimentally available multipartite en- 
tangled states is witnessing a steady growth, stimulated by 
remarkable achievements such as the generation of eight- 
qubit W states GL 01, ten-qubit GHZ-like states fl 0] and 
various sizes of cluster states H]. Very successful lin- 
ear optics experiments have been conducted exploring the 
many entangled classes of up to four-photon quantum cor- 
related states [Jil Ojl [12] and such possibilities are foresee- 
able in other physical systems as well H 1 3fT . Very recently, 
this range of interesting states has been enriched with the 
experimental generation of six-qubit symmetric Dicke states 
and their evaluation in multiparty quantum networking proto- 
cols together with the exploration of a substan- 
tial part of the Dicke class of entangled states. Dicke states of 
n qubits are defined as the eigenstates of both the total spin 
operator J 2 = J2k=x y z^k anc * its ^-component J z , where 

Jk = (1/2) Y^j=i &l an d {&ii&yi & J z} is m e set of Pauli ma- 
trices of qubit j. The model from which this class of states 
arise has been the focus of extensive investigations regarding 
super-radiance, quantum phase transition, correlation and en- 
tropic properties ifTvtl - A Dicke state of n qubits and k excita- 



tions \Dn ; ) is given by 

\D^) = - 
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where is the binomial coefficient and Pi is the set of all 
distinct permutations of 0's and l's. Symmetric Dicke states 
having k = n/2 are particularly interesting in virtue of the 
fact that they are associated with the largest eigenvalues of the 
set of observables { J 2 , J z }. Important theoretical studies fllll 
and experimental progress in characterizing these states have 
been accompanied by the development and optimization of 
special tools explicitly designed in order to "detect" the pres- 
ence of genuine multipartite entanglement (GME) in a state. 
We define GME as states in which all subsystems are entan- 
gled with each other. In this context, a remarkable contribu- 
tion has come from the introduction and use of entanglement 
witness operators (see [kl LL9D and references within), which 



with often only modest experimental effort allow for the dis- 
crimination between separable, biseparable and fully GME 
states without requiring the complete knowledge of the state at 
hand. In fact, it is usually the case that the complete detection 
and quantification of entanglement in a state requires knowl- 
edge of the full density matrix of the system. As this may not 
always be possible, entanglement witnesses provide us with 
viable ways to detect entanglement through only partial in- 
formation. Various forms of witnesses have been formulated 
recently, the most prominent being based on the use of state 
fidelity lEoll and collective spin-qubit operators i2lH22ll . 

It is very important to study the resilience of such char- 
acterization tools to the influences of unavoidable interac- 
tions between the constituents of a given system and their 
surrounding environment. These give rise to phenomena of 
dissipation and decoherence which have a negative impact on 
the entanglement content of a state. Knowing beforehand 
how a chosen method for GME-characterization is able to 
cope with such spoiling mechanisms is not only interesting, 
but also pragmatically useful. It allows one to make predic- 
tions about the performance of a setup and to determine in 
which direction technological progress should be made, in or- 
der to circumvent noise and reliably reveal quantum effects 
for fundamental studies and applications in quantum infor- 
mation tasks. Our work is performed precisely in this im- 
portant direction. We concentrate on the class of symmetric 
Dicke states of n qubits, which exhibit interesting quantum 
properties and are usable resources for quantum networking 
tasks such as quantum secret sharing, telecloning and open- 
destination teleportation OH [3] . We study the behavior of a 
variety of methods for revealing Dicke-class GME, including 
fidelity- and collective spin-based entanglement witnesses, as 
well as less-explored but valuable tools. We provide many 
setup-independent results that can be adapted to those exper- 
imental situations where local measurement settings can be 
reliably and easily arrang ed (as in Refs. HI E S). This 
is a crucial requirement of the detection schemes addressed 
throughout this work. We also discuss feasible techniques for 
the improvement of the resilience of GME-detection to the in- 
fluences of noise. 

The remainder of the paper is organized as follows. Sec- 
tion U briefly introduces the types of noisy channels consid- 
ered in this study. Section [II] studies the behavior, under noisy 
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mechanisms of collective- spin based entanglement witnesses 
of the form experimentally implemented in Refs. |l5|,LL6] and 
proposes an original way to gain robustness. In Section |III1 
we extend our investigation to more common fidelity-based 
entanglement witnesses. These are subjected to filtering oper- 
ations in a way that stretches their tolerance to environmental 
effects and allows for a faithful detection of GME. We ap- 
ply our techniques also to the class of W states and highlight 
some unexpected differences among the channels. In Sec- 
tion |IVJ the GME properties of symmetric Dicke states are 
highlighted by exploiting the interesting network of entangle- 
ment shared by reduced two-qubit states obtained by tracing 
out n — 2 elements. SectionlVlapproaches the problem of reli- 
ably discerning the entanglement class of a given noisy state. 
In Section|VIlwe explore ways to reveal the behavior of multi- 
qubit quantum coherence under the influence of noise, while 
Section PvTO briefly summarizes our findings. Finally, some 
technical details are presented in two appendices. 



I. DECOHERENCE MODELS 

In our study, we address a selection of physically rele- 
vant multi-qubit noisy channels affecting the class of states 
in question. Our choice encompasses a wide range of possi- 
ble mechanisms that are likely to affect a given experimen- 
tal setting designed to achieve multipartite entanglement. We 
use an effective picture for the action of a completely pos- 
itive trace-preserving map given by the operator- sum repre- 
sentation [20]. Within this formalism, a single-qubit noisy 
process is described by a set of Kraus operators {K^}, satis- 
fying the completeness property ^ K^K^ = i, such that, 
calling po the initial state of a qubit, its evolution is given by 
Pch = $ch(po) = KnPoKl> In what follows, for the sake 
of convenience, 7 indicates the characteristic channel's influ- 
ence rate, regardless of its specific nature and ch is a label for 
the channel. 

We start by considering the Kraus decomposition of a zero- 
temperature amplitude damping (AD) mechanism 

K d = |0> <0| + e-^ 2 |l) <1| , Kf = Vl-e-7|0) <1| . 

(2) 

This physically corresponds to an energy dissipation process: 
the system undergoing AD has a finite probability e -7 to lose 
an excitation (here, 7 is an effective dimensionless rate char- 
acterizing the whole process). 

The second process we consider is represented by pure 
phase damping (PD) (or dephasing), which is a phase- 
scrambling and energy-preserving mechanism described by 
the two operators 

*e=jnpi, gf.jEp,,. (3 ) 

It is easy to see that the action of PD on a single-qubit density 
matrix is to exponentially decrease the off-diagonal terms (at 
an effective rate 7), leaving the populations unaffected. Fi- 
nally, we consider a depolarizing channel (DP), which (with 



probability 7) mixes a given one-qubit state with the maxi- 
mally mixed state 1/2. Its action is given by the four-operator 
Kraus representation 

K d P = K? = ^ (k = 1, 2, 3), (4) 

which gives pd p = (l— / y)po-\-jt/2. Therefore, the effect of a 
DP channel is to effectively add white noise to a given single- 
qubit state. This correspondence will become useful in our 
study. The above situation is different from that of a collective 
DP mechanism, which would add white noise to a multipartite 
state of n qubits p mp , leading to (1 — j)p mp + jt/2 n . 

When considering n-partite registers, the effects of equal 
noisy channels, each affecting an individual qubit of the sys- 
tem, can be accounted for by considering q n n-party tensor 
products of Kraus operators /C^'s, where q is the number of 
channel operators of a single-qubit channel. That is 

q n 

Pch = ^£fpm P Kf\ (5) 
j=l 

Our task is now to provide an analysis as general as possible 
of the effects of such environmental channels on a variety of 
experimentally viable tools for multipartite entanglement. 

II. COLLECTIVE-SPIN BASED ENTANGLEMENT 
WITNESS 

Collective- spin operators are useful tools for the inves- 
tigation of GME. Spin- squeezing inequalities fall into this 
class 1 21] and have been extensively studied as well as exper- 
imentally implemented. More recently, collective- spin opera- 
tors that are not directly related to spin squeezing have been 
formulated and shown to be particularly effective when sym- 
metric, permutation invariant states are studied. One can con- 
struct the witness operator i25h 

W n = b bs i-{J 2 X + J 2 y ), (6) 

where is the maximum expectation value of S n = + Jy 
over the class of biseparable states of n qubits. Finding 
(Wn) < for a given state implies GME. The biseparable 
bound, bb s can be numerically calculated (see Appendix A 
and Ref. 1 25]). The witness can be implemented with only 
two local measurement settings making it experimentally ap- 
pealing and realizable in many physical settings (linear optics, 
circuit or cavity quantum-electodynamics). In particular, this 
tool has been used in lUlLUS for the case of four and six-qubit 
states. 

Quite often, Eq. © fails to detect GME in non-ideal sym- 
metric Dicke states which have been affected by noise at their 
generation stage. Through a suitable modification, as exper- 
imentally demonstrated in 1151 . we can provide such a wit- 
ness with greater flexibility in detecting GME in noise af- 
fected symmetric Dicke states. Let us introduce the gener- 
alized collective- spin witness 

W s n (a) = Ma) - S n (a) (qgR), (7) 
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FIG. 1: (Color Online) (a) Collective- spin based entanglement wit- 
ness (S n (a)) for an AD-affected \Df ] ) state. We show the lines FIG 2 : (Color Online) (a) (S n (ot)) for a DP-affected \D^) state. 



corresponding to 7 G [0, 0.3] and highlight the cases where the mod- 
ified entanglement witness turns out to be advantageous: each full 
line is such that (WJ (0)) > with ( WJ (a < 0)) < in a region of 
values of a. The behavior of the biseparability bound is also shown, 
(b) Same as panel (a) but for an AD-affected \D^) state. We show 
the lines corresponding to 7 G [0, 0.16]. In both panels, the shading 
highlights the biseparability region. 



We show the lines corresponding to 7 G [0,0.16] and highlight 
the cases where the modified entanglement witness turns out to 
be advantageous: each full line is such that (VV|(0)) > with 
(W|(a < 0)) < in a region of values of a. The behavior of the 
biseparability bound is also shown, (b) Same as in (a) but for n = 6. 
In both panels, the shading highlights the biseparability region. 



where S n (a) = J 2 + J 2 + aJ 2 . Here, we shall discuss how 
Eq. (|7]) offers more robustness to the noisy channels intro- 
duced in Sec. U than the standard Eq. ©. We notice that the 
bi-separability bound is now a function of new parameter a. 
Using numerics it is seen that in general b bs (a) < b bs (0) for 
a < 0, which implies that at a ^ the threshold for detection 
of GME is lowered. Consequently, we restrict our study to the 
case of negative a. 

Let us start with an AD channel and its effect on the state 
\D^^ 2 ^). We call p a d the channel-affected version of this 

state, p a d = $ a d(|^i n ^)(£>i n ^ |), calculated as described 
in Sec.U One finds 



Tr[S n (a) Pad ] = 'i + ^(n- 1)(1 - e^f 



4 4 v 



(8) 
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This formula is easily proven by noticing that, under a given 
noisy channel, the set of Pauli matrices of qubit j changes as 
°l ch = E u Kl&iK u . For AD, this leads to 0] 



a k,ad — e G 



k ' 



a 3 - 



(i 



(9) 



where k = x,y. From these, it is easy to check that Tr[<r®y <g> 

l® (n - 2) p arf ] =ne"V(2n-2) while Tr[«7f®i 8( "" 2 Va<i] = 
( 1 — e 7 ) 2 — e 7 / (n — 1 ) . By using the explicit decomposition 

of J| given in Appendix A it is straightforward to see that 



Tr[(Jj + J, 2 Kd] = I (2 + ne- 7 ), 

KlJzPad] = \ [n(n - 1)(1 - e- 7 ) 2 + n(l - e~ 2 ^ 



(10) 



which leads directly to Eq. ([5]). We can infer two conse- 
quences. From the first identity of Eq. (ITOl) we see that for 
7 > \n[n 2 /(4b bs - 2n)], we have b bs (0) > Tr[(J 2 + J 2 )p ad }. 
This signals the failure of Eq. © and quantifies the amount 
of AD noise a given experimental set-up can tolerate. The 



second identity of Eq. (ITOl) shows that, different to a pure 
symmetric Dicke state, (J 2 ) ^ for p a d and strongly de- 
pends on 7. From Eq. ([5]) one finds that the dependence of 
the expectation value of S n (a) on a is linear, with a gradi- 
ent coefficient determined by (J 2 ), justifying the inclusion of 
such a term in the witness. For a decohered state such that 
frbs(O) > Tr[(J 2 + Jy)pad]> the inclusion of the z-dependent 
term may be able to pull the expectation value of VV^(a) be- 
low zero, for a set value of a. This is possible if (J 2 ) is 
smaller than the gradient coefficient of the tangent to b bs (a) 
at that point. In Fig. [T] we show two instances of such a pos- 
sibility for the cases of n = 4 and 6, for which bb s (0) ~ 5.23 
and 11.018 respectively (see Appendix A and Refs. fl5[ [l6ll). 
Let us discuss the n = 6 case: At 7 > 0.116 the witness in 
Eq. © cannot detect GME in any AD-affected |L>f } ). How- 
ever, we find instances of decohered states that can be detected 
as GME via VV|(a < 0), although 7 > 0.116. Such a pos- 
sibility goes beyond the example given here. In fact, we have 
strong numerical evidence that for a = 0, GME cannot be re- 
vealed for 7 > 0.076 for an AD-affected | ) . However, for 
a G [—5,-0.1] such states are detected as GME via our mod- 
ified witness VV| (a) despite the fact that they correspond to 
7 G [0.18, 0.20]. Although a proof for any n is difficult (due 
to the computational problem of quantifying b bs (a) for large 
registers of qubits), this shows robustness of VV^ (a) under the 
influences of AD channels. 

A similar analysis conducted with respect to DP noise leads 
to the relation 



'k,dp 



(11) 



which tells us that the (a® 2 



Jn/2) 



n-2) 



) are given by the 

expressions valid for a pure \D)a^ 2 ^) state, multiplied by 
(1 — 7) 2 . Following the same procedure as in the AD case 
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leads us to 



Tr[(4 2 + ^p] = ^ + -(l-7) 



Tr[J^ P ] = -7(2-7), 



(12) 



from which the expectation value of the collective- spin based 
entanglement witness can be determined. In this case too, the 
effectiveness of using VV*(a) is clearly revealed by the ex- 
istence of situations where GME at non-zero a is detected, 
as shown in Fig. O while this is not the case if the standard 
witness is used. 

Finally, we address the PD channel, for which we have (22 



J k,pd 



°k (k = x,y,z), 



(13) 



where we have introduced the Kronecker-delta function Sk z 
which is equal to 1 for k = z and otherwise. This clearly 
implies that a symmetric Dicke state affected by PD noise will 
keep the pure-state property Tv[J^p p( i] = 0, for any n and 7, 
therefore making the strategy highlighted so far ineffective. 
The expectation value of the witness operator reads, regard- 
less of a, as 



Tr[S n (a)p pd ] = - 



-27 



(14) 



which is unable to detect GME as soon as 7 > 
\ ln[n 2 / (46fo s — 2n)]. It remains to be checked whether, for 
instance, an approach such as the one used in Ref. Hfl 
could be used in order to build up a resilient entanglement 
witness in this case as well. 



III. FIDELITY-BASED ENTANGLEMENT WITNESS AND 
FILTERING 

A. Symmetric Dicke States 

We now investigate fidelity-based witnesses 0, [l9] . De- 
spite generally requiring more local measurement settings 
than collective- spin based witnesses, they offer a more spe- 
cific state characterization. Assuming that the state to study 
remains close to the Dicke class, the fidelity-based entangle- 
ment witness with the following general form can be used to 
detect GME 



W n = cJ-\D^)(D^l 



(15) 



where c n is the maximum overlap between \D ( ^ L ^ 2 - > ) and any 
possible biseparable state of n qubits. Quantitatively, c n = 
n/ (2n — 2) for the class of symmetric Dicke states |[l9l l25ll. 

When the AD-affected symmetric Dicke state is studied, 
based on our previous considerations, it is straightforward to 
find that 



Tr[W nPad 



n 



2n 



(16) 



regardless of n. For GME detection we require Tv\W n p a d\ < 



0, which is guaranteed for 7 < (2/n) m[(l/c n )]. For n — 4 
(n = 6) the witness is implementable with 9 (21) local 
measurement settings (see Appendix B and Refs. [0 [l6ll ). 
The limiting amount of sustainable AD noise for n = 4 is 
7 = 0.203 and decreases as larger Dicke states are examined 
(we have 7 = 0.170 for n = 6, see Fig. [3] (a)). Fig. @](a) 
shows this "pulling-back" effect for n = 4, .., 50. As previ- 
ously done, we seek to devise experimentally-friendly tech- 
niques to detect GME for larger amounts of noise. To achieve 
this task, we apply filtering operations GHHH] to the fidelity- 
based entanglement witness. The local nature of the filters 
cannot alter the amount of entanglement in a state. However, 
they may allow for an increase in the noise allowed before a 
witness starts failing to detect GME. For n-qubit states, we 
thus construct filters of the form 



(17) 



where JF^'s are local invertible operators. This technique has 
already been investigated for a variety of states in Refs. (US 
[sl H3] and here we extend and generalize its use to arbitrarily- 
sized symmetric Dicke states. 

We must be careful to ensure that the correct normal- 
ization is taken. To this end we impose the constraint 
TrfWnl^TrfW^ 7 ] with me filtered witness operator. This 
is sufficient to ensure that the expectation value arising from 
the new witness is comparable to the expectation value from 
the unfiltered one. Therefore 



w = 



Tr[^W n J^t] 



(18) 



Filtering may lead to an increase in the number of required 
local measurement settings. However, there are indications 
that the class of filters given by B24I1 



1 

yj 



(19) 




FIG. 3: (Color Online) (a) Expectation value of an entanglement 
witness for a symmetric six-qubit Dicke state undergoing AD at a 
dimensionless rate 7. The solid curve corresponds to the unfiltered 
witness while the dashed curve is for the filtered one. By filtering 
one can increase the range of 7 where the witness is still able to de- 
tect GME. (b) Expectation value of the entanglement witness for a 
symmetric six-qubit Dicke state undergoing PD type of noise. The 
solid curve is for the unfiltered witness while the dashed curve cor- 
responds to the filtered case. No advantage is achieved, in this case, 
upon filtering. 
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(a) 



(b) 
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FIG. 4: (Color Online) (a) Tr[W n /w] plotted against the dimen- 
sionless AD rate 7 for n = 4 — > 50. Solid (Dashed) lines are for 
unfiltered (filtered) witnesses. As n grows, the filtering-induced gain 
is reduced, (b) Values of 7 at which the filtered witness surpasses 
the threshold value — 10 _3 and the unfiltered one starts to become 
non-negative with n. The lower (upper) curve is for the unfiltered 
(filtered) case. The advantage achieved upon filtering is clear as it 
slowly decays as n increases. 



with i/j a positive real number, guarantees for the small- 
est possible number of necessary measurement settings. A 
direct calculation, performed using (for convenience) yj = 
y (Vj = 1, .., n), leads to 



Tr[W£p ad } = 



(2 n n-2n + 2)e- 



n (y 2 



(2n - 2)y n 



L 2n - 2 



(y 2 



1) 



n/2 



(20) 



This is then minimized with respect to the filtering parameter 
y for any set value of 7, ensuring that the expectation value 
of the witness would not depend on y at all B24I1 . The effect 
of this procedure can be clearly seen in Fig. 0(a): the local 
filtering increases the amount of decoherence tolerated while 
the witness is still able to detect GME in a symmetric Dicke 
state <28h. 



The actual detection of GME via negativity of a witness 
operator clearly depends on the errors associated with the ex- 
perimentally determined value of Tr[VV^p a( J: a small (albeit 
negative) value is likely to be covered up by the corresponding 
error bar. Although the quantification of any acceptable lower 
bound for significant GME detection is a setup-dependent is- 
sue, based on current linear optics implementation, it is rea- 
sonable to expect that, upon collection of a sufficiently large 
sample of data, values of Tr[W^ P ad] 00 — 10 can still be 
discerned from zero. We thus fix a threshold of this order 
of magnitude and seek the smallest value of 7 at which the 
filtered witness (for a given n) surpasses it. This provides a 
practical lower limit to the mathematically rigorous perfor- 
mance of such tool. Fig.|4](b) reveals the advantage acquired 
upon filtering for n = 4 — > 20: the gap with respect to the 
unfiltered case is quite considerable and marks the success of 
this strategy. Quantitatively, the upper points in Fig. 01 (b) are 
well fitted by the function ln[( I ^)i a^] with ai = 1.54, 
b x = -1.54 and ci = 4.73. 

We now study the detection of GME via a fidelity-based 
witness in symmetric Dicke states affected by PD type of 



noise and find 



Tr[W nPpd ] 



2n 



/2 



n/2 

£(^/ 2 ) 2 

k 



-7(n-2fe) 



(21) 



One can easily find that the general qualitative features high- 
lighted for the case of the AD-related study of an unfiltered 
witness hold in this case as well: as n grows, the values of 
7 at which the expectation value of Eq. (I2T1) become posi- 
tive are quickly pushed towards zero. However, an impor- 
tant remark is due: different to the case of AD noise, for 
PD the filtering technique employed above does not increase 
the range of tolerated noise. In fact, one can easily see that 
the net effect of the application of filtering operators is that 
Tr[W^p p d] = y n TY[W n p p d\. Therefore one cannot shift the 
value of 7 at which GME is unambiguously revealed. This is 
shown by the dashed line in Fig.[3](b), where we see that the 
expectation value of the filtered operator tracks its unfiltered 
version in the negative semi-space |29|]. The reason behind 
such a behavior is clearly understood by the noticing that, 

n^nP P d} = Tr[W n ^^^f |^/ 2 ))<^/ 2 )|^ dt ^, 
3=0 

(22) 

where, as defined in Sec.H K> P j 's are the n-qubit Kraus oper- 
ators for the whole PD-affected register. We can thus interpret 
this equation as the expectation value of the unfiltered witness 
over a symmetric Dicke state affected by new Kraus opera- 
tors, each given by ^K pd . It is matter of direct calculation to 
see that the resulting density matrix p p d is simply y n p p d, thus 
demonstrating our claim. This is obviously not the case for 
an AD channel, in virtue of the different form of single-qubit 
Kraus operators required in that case. In order to exclude any 
limitations induced by the choice of the specific form of filter- 
ing operator, we repeated our analysis using general invertible 
operators, still finding no improvement upon filtering. 

A very similar situation is encountered for the case of DP 
channels acting on the qubits of the register. In this case, the 
number of Kraus operators involved in the evolution of a given 
state grows as 4 n . This makes any approach to the problem 
intractable for large n. Therefore, we have not been able to 
produce a general formula for the fidelity-based entanglement 
witness nor a simple argument to explain why, in this case as 
well, filtering is non-effective in providing noise-robust GME 
detection. In Figs. \5\ (a) and (b) we give evidence of this for 
n = 4 and 6, comparing the unfiltered witness with a few 
instances of filtering. Moreover, panel (c) shows that the fil- 
tered witness has an absolute minimum at y = 1 in its region 
of negativity, demonstrating the ineffectiveness of filtering. 



B. W-States 

With minimal changes to the analysis performed with re- 
spect to the symmetric Dicke states, we can also study LD^ ), 
which are commonly referred to as n-qubit W states 1 1], thus 
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(a) 



(b) 



Tr[W iPdp ] 
0.6 




Tr[)V 6Wp ] 

0.6r 





y 5^ 



IV. REDUCED-STATE BASED ENTANGLEMENT 
WITNESS 



We now consider a different manifestation of GME in sym- 
metric Dicke states based on the observation of entanglement 
residing in the two-qubit reduced states of \D ( f h ^ 2 ^). Tracing 
out n — 2 qubits, we find 

q = a n \^){^\ + ii^^[|00)(00| + |11)(11|] (26) 



with a n = n/[2(n - 1)] for n > 4 UM and |^ ± ) 
(|01) ± \10))/V2. Here, g has fidelity HI (F^+) 
Tr[\ip+)(ip+\o\ = a n with respect to \^ + ). Using the fidelity- 
based entanglement witness 



^spe( 
L19] 



1 



(27) 



FIG. 5: (Color online) (a) Expectation values of filtered and unfil- 
tered fidelity-based entanglement witnesses for DP-affected l^ 2 ^) 
plotted against the characteristic DP rate 7 (dimensionless). Each 
dashed line is for a filtered witness with y £ [0, 3] at steps of 0.1. 
(b) Same as in panel (a) but for \D^). (c) Filtered entanglement 
witness for DP affected \D^) against both y and 7. In its negativity 
region, the witness is minimized at y = 1. 

showing the versatility of the techniques employed in this pa- 
per. We consider the fidelity-based witness IU9I1 

n - 1 



-1 



n 



\dP)(dP\, 



(23) 



where in order to distinguish this case for the previously 
treated one, u c h is used to indicate the density matrix result- 
ing from the application of channel ch to a pure n-qubit W 
state. Under AD, the expectation values of Eq. (f23l) for the 
unfiltered and filtered cases are given by 



Tr[W n uu ad ] 



11 -e~ 7 , 

n 

[2 n (n- 1) -n][(n- 1)(1 



n[(n - 1) (y 2 + if 



ny A 



(24) 



The effect of filtering once again increases the amount of tol- 
erated noise for the witness to still be able to identify GME. 
Considering PD, we see that the addition of filtering has no 
benefit for the witnesses (for the same reasons explained in 
the discussion put forward in Sec. lIII At . We find the expecta- 
tion value of the witness to be 



Tr[W n u; pd } 



(n- l)e 



-27 



(25) 



Under appropriate conditions, these results are in agreement 
with the study reported for n = 4 in |24j]. As for the DP 
channel, our investigation reveals that, at least for n = 4 and 
6, a small advantage is gained by filtering the fidelity-based 
witness, although we do not have a closed analytical form for 
any n. 



we have (W r ) = 1/2- Thus, a fidelity (F^+) > 1/2 

detects the presence of entanglement in the two-qubit reduced 
state q. This is always possible for all choices of pairs of 
qubits in |Z)i n ^), as a n > 1/2 V n. However, clearly, this 

may not be true for other states, including \D ( f 1 ^ 2 ^) subjected 
to the noise channels outlined in the Section U In general, 
we define two qubits of a multipartite state | <j>) as connected 
if their reduced density matrix is such that (F^+) > 1/2. In 
this sense, a symmetric Dicke state gives rise to a connected 
set of reduced states. For any given n, one can construct a 
graph having qubits j = 1, .., n at its vertices. Two vertices 
are joined by an edge if and only if they are connected in the 
sense explained above. According to this definition, symmet- 
ric Dicke states give rise to fully connected (complete) graphs, 
as shown for n = 4 and 6 in Fig. [6] (a). 

In order to relate the witness W r to observables in an ex- 
periment, we decompose the fidelity {F^+ ) into expectation 
values of Pauli operators as follows 



(F 2 j,+ ) = -(l + Tr[a x <g)a x Q\+TY[ay <S> cr y g] - Tr [a z ®a z g}), 

(28) 

showing that only the local measurement settings a® 6 , a® 6 
and af 6 are required for an experimental implementation. In 
fact, any two-qubit correlation Tr [a l k <S> (J 3 k g] can be obtained 
from the corresponding data. This puts the method described 
here on an equal footing with the collective- spin based witness 
described in Sec. Ill in terms of required experimental effort. 

We now investigate how this method copes with noise af- 
fecting the class of symmetric Dicke states. By using Eqs. ©, 
(TUt and (fT3l) . together with Eq. (l28t , it is straightforward to 
show that 



- 2 T(a„ + e^(l+a„)-l), 



(^ + )* = «n(7-l) 2 - i (7-2) 7 , 



(29) 



pd 



-(l + e- 2 ->„. 
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FIG. 6: (a) Fully connected graphs for \D± ) and \Dq ). The ver- 
tices represent qubits and the edges represent the presence of entan- 
glement within their reduced- state, as detected by using (W r ). (b) 
Effects of noise on the detection of entanglement in the reduced two- 
qubit states of \D^^ 2 ^) using (VV r ). Here, the solid (red), dashed 
(blue) and dotted (green) lines correspond to (F^+) for the AD, DP 
and PD channels respectively. The shaded grey area corresponds to 
the region in which (VV r ) fails to detect entanglement and the corre- 
sponding graph in (a) becomes completely disconnected. 

In Fig. [61(b) we show the behavior of (F^+) c/l for the three 
channels. As soon as 7 > ln[(l + a n + (a^ + 6a n -3) 1 / 2 )/2], 
7 > 1 - (4a n - l)- 1 / 2 and 7 > -|ln[(l - a n )/a n ] re- 
spectively, one finds disconnected sets in an AD-, DP- and 
PD-affected \D { n /2) ) state using (W r ). These thresholds are 
larger than those corresponding to the use of a collective- spin 
based entanglement witness for any value of n that we could 
quantitatively consider (see discussion in Sec. (TTJ), as a result 
of the smaller dimension of the states being tested. 



V. STATE DISCRIMINATION VIA CHARACTERISTIC 
OPERATORS 

The variety of ways multipartite entanglement can be 
shared by an n-qubit register requires ways to determine if an 
experimental state belongs to one of the known classes of en- 
tanglement. This can be achieved by relying on the formalism 
of state discriminators 111 111, which have been experimentally 
implemented and used in order to assess the classes of four- 
qubit entangled states 11 ill . Here, we study the reliability of 
such methods for symmetric Dicke states suffering effects of 
noisy channels. 

Consider multi-qubit operators having \D^^ 2 ^) as a non- 
degenerate eigenstate associated with the largest possible 
eigenvalue of the operator's spectrum. We call such oper- 
ators characteristic of the state \D^^ 2>} ). At least one op- 
erator having these features exists, i.e. the fidelity operator 
\D^/ 2 - > )(D ( £ l ^ 2 - > |. However, this is not the only option and 
other characteristic operators can be designed, requiring far 
less local measurement settings than the fidelity one. A sys- 
tematic approach would require the decomposition of the fi- 
delity operator into tensor products of single-qubit Pauli oper- 
ators (see Appendix B). Out of them, only the genuine n-qubit 
correlators should be selected to be combined together in a 
way so as to construct a proper characteristic operator. For 
instance, for n = 6 qubits prepared in \D^) we can build 



up the Bell-Mermin operator JH B D (3) = T,k=x, y ^V 20 
with 6 k = a\ [af 5 - Ei Pi(°f 3 ® °T - &k ® 

This is characteristic for \D^), which is an eigenstate with 

associated eigenvalue 1 (the maximum within 23 ( 3 )'s spec- 

6 

trum). However, if we simply take the negative terms V (3) oc 

6 

- T,k=x, y & l ® Pi{?t 3 ® & f 2 ) in this expression, we find 
that \D^) is still an eigenstate of maximum eigenvalue (en- 
forced to be 1 upon renormalization). Representatives of any 
other class of entanglement will achieve expectation values 
smaller than 1. This can be used for effective entanglement- 
class discrimination. For instance, six-qubit GHZ states trans- 
formed by local unitary operations (LU) or stochastic local op- 
erations supported by classical communication (SLOCC) Il20ll 
yield expectation values no-larger than 0.833. Thus, if an ex- 
perimental state p exp of n = 6, thought to be close to the 
symmetric Dicke family, gives {T> Pexp ) > 0.833, one can ex- 
clude any GHZ-like character. This can be adapted to any 
other class of GME states. Here, without affecting the gen- 
erality of our study, we concentrate on the discrimination be- 
tween noise-affected symmetric Dicke states and the n-qubit 
GHZ class. We thus construct the streamlined characteristic 
operators 

p D <„ /2) =- n £ *i® y, w~ 3 ® <^ 2 )> oo) 

k=x,y I 

where M is a normalisation factor taken so that the eigen- 
value corresponding to \D^n^ 2 ^) is 1 and found by notic- 

ing that -<Z)i n/2 Vr~ 2 ® *f 2 \D { n /2) ) = n/(2n - 2) for 
k = x,y and that the number of possible permutations in- 
volved in Eq. is C^ 3 , so that M = 2/[n(n - 2)]. We 
thus see the effects that a channel has on the expectation value 
of discrimination operators. For AD, PD and DP noise, re- 
spectively, are 

Tr[P D (. /2) p ad ] = -e- I ?(e^-2), 

Tr[D D(nn/2)Ppd }=e-^\ (31) 

Tr[P D (. /2 )Pdp] = (l-7) n . 

The first two equations can be understood by using arguments 
analogous to those valid for collective- spin witness opera- 
tors. The third identity of Eq. (I3T1) is clarified considering 
the analogy between depolarizing channel and single-qubit 
white noise. In Fig. [7] we show the behavior of Eqs. (I3T1) 
for n = 4, 6, 8 and 10 as the amount of the respective noise- 
influence increases. We also show the bound associated with 
the n-qubit GHZ class. Thus, the shading in Fig. [7] represents 
the region where discrimination between one of the channel- 
affected Dicke states and the GHZ class is not possible. After 
a certain noise strength, we can no longer determine if the 
state being studied is a noise-affected Dicke state or a state 
from the GHZ class. The decay increases more sharply with 
larger n and the DP channel has the worst effects. We are 
working on the design of a strategy based on filtering opera- 
tions which might allow one to gain robustness of this discrim- 
ination procedure against noise. However, the potential of this 



FIG. 7: (Color online) State discrimination via characteristic oper- 
ators for symmetric Dicke states with n — 4, 6, 8, 10. The contin- 
uous lines show the behavior of AD-affected \D^ l ^ 2 - > ), the dashed 
lines are for a PD channel, while the lowest dotted lines are for a DP 
mechanism. Each channel is assumed to be characterized by a rate 
7. The horizontal lines show the maximum expectation values of 
the characteristic operators T> D ( n /2) over SLOCC-equivalent GHZ 
states of n qubits. The shading highlights the regions where discrim- 
ination between symmetric Dicke states and GHZ class is no longer 
possible. 

tool is already seen at the level of practical implementability. 
For instance, the n = 6 version of Eq. (f30b can be decom- 

posed as T> D m =E fe=x ,/M| ®5= 2 (^H ®i= 2 (^- 

H)-TS ®U& + 2a b+T2 ®'= 2 (^-2«Ti) + 5(8) J 6 =2 a{], 
which only requires 10 local measurement settings for its im- 
plementation. Finding out the explicit decomposition of dis- 
crimination operators for general n is, however, a daunting 
problem. 

VI. CORRELATION FUNCTION 

In this Section we shift the focus of our discussion from 
GME to the quantum-coherence properties of the class of 
states under investigation. Our approach here considers the 
expectation value of the n-qubit correlation operator 

C(tf) = (costf<7 fe + sintf<7 j )® n (32) 

with k 7^ j — x,y,z. Through this, one probes the coher- 
ence of each element of the register along a direction (in the 
single-qubit Bloch sphere) lying in the plane formed by the 
unit vectors k and j. The use of such a multi-qubit correlator 
is common in the assessment of the properties of GHZ states 
for quantum metrology purposes |32]: as a result of n-qubit 
coherence, when k = x and j = y the expectation value of 
Eq. (f32l) oscillates with $ at a frequency that depends on n. 
Here, we shall investigate collective coherence in an n-qubit 
symmetric Dicke state by means of Eq. d32b . This requires the 
implementation of a single measurement setting per value of 
$ and is thus an experimentally favorable tool for multipartite 
state characterization. For the sake of definiteness, here we 
concentrate on k = x and j = z, although any other choice is 
equally suitable. 

For a pure |D n n ^ ) state, the expectation value of the multi- 
point correlator is easily found using basic combinatorial ar- 



FIG. 8: (Color online) Behavior of the expectation value of the n- 
qubit correlation function calculated over symmetric Dicke states 
|L>n n/ ^) with n — 4, 6, 8 and 10. A parity-dependent effect related 
to the number n/2 of excitations in the state being studied is seen: 
states with an even (odd) number of excitations give rise to a positive 
(negative) expectation value at = ty/2. Moreover, the position of 
secondary maxima and minima is n-dependent. The beating is an 
effect of quantum coherence in the state. 

guments and the symmetries in the class of states at hand. We 
have 

n/2 

(D("/ 2 )|CW|D("/ 2 )) = ^(-l) fc (C„ fe /2 ) 2 ( C o S ^)"- 2fc ( S in^) 

k=0 

(33) 

In general, Eq. (l33t exhibits an n-dependent oscillatory be- 
havior whose features strongly depend also on the parity of 
n/2, as shown in Fig. [8] The expression corresponding to 
n = 6 has been used in Ref. 1 15] in order to contribute to the 
characterization of | ) . 

Here, we shall study the behavior of (C(#)) when a sym- 
metric Dicke state is subjected to noise effects. By using 
again the expressions of the single-qubit Pauli operators trans- 
formed upon the action of a given channel (see Sec. [II]), one 
can prove that Eq. (l33t remains almost invariant under en- 
vironmental action. The only modification is at the level of 
oscillation amplitudes, which are changed by a 7-dependent 
factor. Explicitly, we have found the following universal form 

n/2 

Tr[C»p c ,]=^(^ /2 ) 2 r^( 7 )(cost?)"- 2fe (sint?) 2fc 

k=0 

(34) 

with ch = {ad,pd, dp} and 1^(7) = e -(t- fe h(i _ 2e~^) k , 

r&(7) = (-l) fe e-(- 2 ^ and Y% (7) = (-l) fe (l - lY 
being the channel- specific factors responsible for the loss of 
coherence in the state. Our claim is that the beating ef- 
fect responsible for the rich oscillatory structures shown in 
Fig. [8] arises only in virtue of the quantum coherences within 
|D n n/2) ). In fact, while (7) - as 7 grows, so as 
to progressively kill any oscillations, r^ d n / 2 (7) ~* 1 anc * 

ri/2 W - (- 1 )" 72 ' with KTh) - Vfc < n/2. This 
implies that, as the coherences in the n-qubit state disappear 
under an AD or PD channel, the modulus of corresponding 
correlation functions simply becomes |(sin$) n |. Thus, al- 
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(a) (b) (c) 




(d) (e) 




FIG. 9: (Color online) We show Tr[C(tf)/w] [panel (a)], 
TY[C($)p pd ] [panel (b)] and Tr[C($)p dp } [panel (c)] against & and 
the respective 7 for a noise affected \D^) state. Parity effects as- 
sociated with the various channels considered are shown. Panels (d) 
and (e) show the AD and PD case, respectively, when the classical 
asymptotic behavior is subtracted. 



though an oscillatory behavior is still kept in these cases, no 
beating is found. A difference can be found between the 
trends corresponding to AD and PD channels: a PD chan- 
nel would progressively destroy the off-diagonal elements of 
a density matrix without affecting its populations. Asymptot- 
ically, this results in a diagonal state whose only non-zero en- 
tries are those corresponding to states having n/2 |0)'s. If we 
now take d = n/2 in the n-qubit correlation operator, we eas- 
ily understand that its expectation value over the PD-affected 
state can only be ±1, depending on the parity of n/2. On the 
other hand, such a parity-effect is absent in the AD case. In- 
deed, asymptotically, this channel would reduce any symmet- 
ric Dicke state to its collective ground state, which can only 
give (C(7r/2)) = 1, Vn. These features are well illustrated in 
Figs. [9] (a), (b) and (c) for the case of a symmetric six-qubit 
Dicke state. 

Such an asymptotic classical behavior tends to mask the 
trend that the beating follows in the changes undergone by the 
state. We have thus stripped the correlation functions from 
such contributions by subtracting their Fourier- series expan- 
sions from the analogous one of Tr[C(#)p a d jP d]. The results, 
which highlight the sole decrease in visibility of the fringe 
of beating induced by quantum correlations, are shown in 
Figs. 0(d) and (e), which show an evident exponential decay 
against 7. Such a procedure is not necessary for the DP, which 
smoothly flattens the correlation function to zero. The iden- 
tification of a general trend against the number of qubits in a 
state is a task made difficult by the ^-dependence of such fig- 
ures of merit: states affected by different channels give rise to 
maxima and minima of Tr [C (#) p c h] located at different values 
of the angle One can extract useful indicative information 
by looking, for instance, at the correlations corresponding to 
d = 0, for a set channel and increasing number of qubits, 
to find that the exponential decay of correlations induced by 
a growing 7 becomes faster for larger n. Overall, the ana- 
lytic expressions provided here for this set of relevant noise 



channels embody a valuable tool for the experimental charac- 
terization of symmetric Dicke states. The theoretical curves 
can indeed be used to fit the points acquired, experimentally, 
by tuning the parameters of a setup in a way so as to properly 
select the direction $ along which one would like to probe 
quantum coherence lfl5h . 

VII. CONCLUSIONS 

We have studied both GME-detecting and state- 
characterizing tools for the class of symmetric Dicke 
states of an arbitrary number of qubits, under the influences 
of noise. Our investigation starts from standard instruments 
such as fidelity-based entanglement witnesses and collective- 
spin operators and unveils some experimentally-friendly ways 
to improve their resilience to noise. Besides its pragmatic 
features, which make our study explicitly devoted to the pre- 
diction and evaluation of the performances of an experimental 
setup, we have revealed interesting characteristics of noisy 
channels and their influence on GME witnesses. For instance, 
we have clearly shown that resilience to PD channels cannot 
be achieved by means of simple filtering operations per- 
formed on fidelity-based entanglement witnesses. In light of 
the steady technological progress in a variety of experimental 
settings, ranging from linear optics to solid-state systems, we 
expect our results to be useful for the purposes of a complete 
characterization of multipartite entangled states generated or 
evolving in the presence of environmental effects. Our results 
should also help in the context of reliably determining the 
quality of experimental symmetric Dickes states for use in 
quantum networking tasks such as quantum secret sharing 
and open-destination teleportation. 
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APPENDIX A 

In this Appendix, we outline the steps involved in determin- 
ing the biseparability bound used in the entanglement witness 
based on collective- spin operators. We start by considering 
the collective operator 

S n = Jl + Jl+aJ 2 z (A-l) 

with aGl. Upon expansion of each collective- spin operator 
and using the "one-vs-(n — 1)" qubit bipartition, we get the 
expression 

^ rjr^ TIOL 1 A /v 

Sn = (- + — n + -(<J 1 x Q x +a 1 y Q y +aa 1 z Q z + R a ), (A-2) 
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where Q k = YJj=2 °k and R a = R x + R y + olR z with R k = 
YTjZl °k E/=j+i = x > V-> z )- We take me expectation 

value of S n over the biseparable state of qubit 1 and system 
2 — > n, so that 

0n> = | + ^ + ^[E (^><4) + ^>(4) + (^)] 

n na 1 . n not 1 /A N 7 . . 

= 2 + T + 2 (a>> - 2 + T + 2 M 2 < X i ( a) = M " } 

(A-3) 

with r = z), {A a } the set of eigenvalues of Cl a , x = 

and analogous expressions for y and z. The diagonal- 
ization of Ct a can be performed numerically for a few rele- 
vant cases, therefore allowing the quantification of the upper 
bound for biseparability associated with this splitting. For the 
specific cases of n = 4, 6 and 8, the "one-vs-(n — 1)" splitting 
gives the largest bound among all other possible bipartitions, 
thus providing the threshold bb s (a) for the detection of GME 
(see SecUB). Quantitatively, we have b bs (0) = 5.23, 11.018 
and 18.83 for n = 4, 6 and 8 respectively. 



APPENDIX B 

We discuss the local Pauli decomposition of the fidelity op- 
erator \D^^ 2 - > )(D ( £ l ^ 2 - > | and provide a useful method for re- 
ducing the number of local measurement settings experimen- 
tally required to measure it. To decompose a given n-qubit 
projector \^>)((j>\ into local Pauli form, we use the relation 



!</>> (</>! = ^ Yl ®<7i 2 ® .. ®£i n ), (B-i) 

ii ,i 2 ...i n 

where we have introduced the correlation tensor C^,^...^ = 
& il (8) &i 2 ® .. <g) &i n |0) and i n e {0, x, y, z}, with &o = 
i. Using Eq. (IB- lb , the correlation tensor corresponding to 

(2) (2) 

\D y /')(D y /'\ has 40 non-zero elements 



k=x,y,z 



2(16"® + 1 <7® + 1 CT® — (T® <7 



®2_ /T ®2^ r ®2_ (T «.2 /T ®2 



(B-2) 

where indicates all distinct permutations of the opera- 
tors. However, such a 40-element decomposition can be sig- 
nificantly streamlined by using the relation (valid for z, j = 

x,y,z) 

E *f 2 -f 2 = ^ + -,) 04 + - - 4 - , 

(B-3) 



through which one can rewrite the ^2 term of Eq. (IB-21) 
in terms of the 9 local measurement settings :&f£ iZ9 [(c x ± 

a^)/^] 04 and [(cr y ia^/V^] 04 . Thus, the symmetries 
present in the state's decomposition have been exploited in a 
way so as to reduce the number of measurements. Using the 

(2) 

compacting techniques outlined here for | D\ ; ) , one should be 
able to apply them to arbitrary sized symmetric Dicke states 
\D^ l ^ 2 - > ) to obtain significant reductions in the number of lo- 
cal measurement settings. 

For example, using Eq. (IB-1K one finds that the correla- 
tion tensor corresponding to \D^)(D^\ has 544 non-zero 
elements 



'320 

T ®2^®2 , fl ®2^ 



l^fffffj" + r z <rf - 1TV® 2 - 1® 2 <t® 2 af 2 
_3_ 
320 1 

<7® 4 <7® 2 + 1® 4 <7® 2 + l® 2 <jf 4 + l® 4 (j| 2 + l® 2 <j| 4 ]) 

^[a® 6 +af + l® 6 -a® 6 ]. 



1® 2 < a® 2 ] + ^[<x x ® 2 a® 4 - afaf 4 - a x ® 4 af 



Using Eq. dB-3l >. together with the relations 



(B-4) 



J2 <rfafaf = ±[(«r x + a, + <x z )® 6 + («x x + ff „ - «r z )« 



+ (ffi _ ffw+(Tz) ®6 + (ffit _ ffi/ _ ffz) ®6 ] _^ (ff |2 (T ®4 +(T |4 ff ® : 

7T 

7r k=x,y k=x,y,z 



5>f 2 <rf - of V® 2 ) = ^{(a, + 2a,)® 6 + (a, - 2a,)® 

7T 

- 4[(a; + a,)® 6 + (a; - a,-)® 6 ] - lOaf 6 - 136a® 6 }, 

(B-5) 

one can rewrite the permutations of Eq. (lB-4b in terms of 
the following 21 local measurement settings: c® 6 ^, [(ct x ± 
cj mz )/V2)®\ l(a y ±a z )/V2}® 6 , [(* Xty ±2* z )/y/5]® 6 , [(a z ± 
2<7x iy )/V5\ 06 and [(cr x ± ^±a 2 )/v^] 06 . Such a reduction 
provides a huge advantage for measuring the fidelity of | ^ 3) ) 
in a given experimental setup. 
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